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Abstract 

We study a dyadic branching Brownian motion on the real line with absorption at 0, 
drift /r G M and started from a single particle at position x > 0. When /r is large enough 
so that the process has a positive probability of survival, we consider K{t), the number of 
individuals absorbed at 0 by time t and for s > 0 the functions u}s{x) := We 

show that Wg < (X) if and only of s G [0, Sq] for some Sq > 1 and we study the properties of 
these functions. Furthermore, for s = 0,uj{x) := ujo{x) = P^(A((X)) = 0) is the cumulative 
distribution function of the all time minimum of the branching Brownian motion with drift 
started at 0 without absorption. 

We give three descriptions of the family s G [0, sq] through a single pair of functions, 
as the two extremal solutions of the Kolmogorov-Petrovskii-Piskunov (KPP) traveling wave 
equation on the half-line, through a martingale representation and as an explicit series ex¬ 
pansion. We also obtain a precise result concerning the tail behavior of K{oo). In addition, 
in the regime where K{oo) > 0 almost surely, we show that u{x,t) := = 0) suitably 

centered converges to the KPP critical travelling wave on the whole real line. 


1 Introduction 

Consider a branching Brownian motion in which particles move according to a Brownian motion 
with drift /r G M and split into two particles at rate /3 independently one from another. Call 
A/^aigP the population of all particles at time t and call Xu{t) the position of a given particle 
u G A/'aii(i). When we start with a single particle at position x we write P* for the law of this 
process. 

In a seminal paper, |ld| . Kesten considered the branching Brownian motion with absorption, 
i.e. the model just described with the additional property that particles entering the negative 
half-line (—oo, 0] are immediately absorbed and removed. We write for the set of particles 

alive (not absorbed) in the branching Brownian motion with absorption and K (t) the number of 
particles that have been absorbed up to time t. The system with absorption is said to become 
extinct if > 0 : Muveit) = 0 and to survive otherwise. We let K(oo) := lim^-^oo K{t) G MU{oo}. 

Depending on the value of fj, one has the following behaviours (see Figure]^ 

Regime A: if /U < — \/2^, the drift towards origin is so large that the system goes extinct almost surely. 

K{oo) is finite and non-zero. 
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Regime B: if —-y/2^ < fi < there is a non-zero probability of survival. On survival, there will 
always be particles near 0 and K{oo) = oo almost surely. 

Regime C: if /r > y/^ there is still a non-zero probability of survival, but the system is drifting 
so fast away from 0 that, on survival, min^g_;\/^jj(j) Xu{t) drifts to -|-oo almost surely as 
t —7- cx); K{oo) is thus almost surely finite. Furthermore, there is a non-zero probability 
that K{oo) = 0. 





Regime A ■. /i < —\/2/J Regime B : — \/2/J < /i < \/2/J Regime C : /x > \/2/J 

Figure 1: The three regimes: The hashed region represents the cone in which one expects to 
find particles. 


The behaviour of K(oo) in regime A {^ < —y/^) has been the subject of very active 
research recently, including a conjecture by Aldous which was recently settled by P. Maillard 
|15| (we discuss Maillard’s results bellow), improving earlier results of L. Addario-Berry and N. 
Broutin [Tj and E. Ai'dekon j2]. Surprisingly, relatively little was known concerning the regimes 
B and C. Our main results in the present work concern the study of K(oo) and of certain related 
KPP-type equations. 


1.1 The tail behaviour of K{oo) 

In |15) . Pascal Maillard has shown that in regime A (^u < —y/^) the variable K{oo) has a very 
fat tail. More precisely he shows that, as z ^ oo, there exists two constants c,c' which depend 
on X such that 


F^[K{oo) > z] ~ 


c/{zlog{z)‘^) 


for /X = -\/2)d, 

for /X < —y/^ where a(/x) 




In regime B {—y/^ < /x < y/^) it is clear that K(oo) = oo on survival so one would 
essntially condition on extinction to study the tail behaviour of A'(oo). This is outside the scope 
of the present work. 

In regime C (^ > y/^), however, K{oo) is almost surely finite. We introduce for s > 0 and 

X > 0, 

a;.(x) a;.(0) = s. (1) 

When s G [0,1] this quantity is the generating function of K{oo). We show that Ws(x) is finite 
for some values of s larger than 1: 

Theorem 1. In regime C (fi> y/^), there exists a finite sq > 1 depending only on n/y/fi such 
that 

1. For s < So, Ws(x) is finite for all x > 0, 

2. For s > So, Ws(x) is infinite for all x > 0. 
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The functions x i— )• 0 Js{x) are increasing for any s G [0,1) and decreasing for any s G (l,So], 
converging to 1 when x —)> oo, and one has w^q(O) = 0. 

Furthermore, one has for n large 


r^[K{oo) 


2s'^ni-^TTl5{sQ - 1 ) 


( 2 ) 


Using the branching strnctnre and a simple conpling allows to relate the ujs{x) with each 
others, as in the following resnlt (which is already present in |17[ I15]l: 

Theorem 2. In regime C (gL> \/2^), 

1. For each s ^ [0,1), ujs{x) = ujo{x + lOq^{s)), 

2. For each s £ {1, sq], uJs{x) = uJsq{x + U!T^^{s)). 

We do not have an explicit expression for sq as a fnnction of yL/\/j3, bnt we can evalnate it 
nnmerically with a good precision as shown in Fignrej^ In the critical case /r = \/2[d, we obtain 
So = 1.3486.... 



t/'/P 

Fignre 2: Nnmerical determination of sq and s^^j[if as fnnctions of 
We can prove the following property: 

Proposition 3. sq is an increasing function of ti/\/]3 and furthermore sq ~ cijf /fi for some 
constant c as [ify/~^ —t oo. 


1.2 Distribution of the all-time minimum in a branching Brownian motion 

The probability that K{oo) = 0 for a system started from x, is also the probability that the 
all-time minimnm of a fnll branching Brownian motion with drift /i started from zero does not 
go below —x: 


ui{x) := UJo{x) = P^[iF((X)) = 0] = P'’[min min Xu{t) > —x], (3) 

t>o ueA/'aiih) 


This quantity, of conrse, is not trivial only in regime C {[i > \/2i3). Then, since 


lim min X^ft) = -|-oo 

i^OO ugATallh) 
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almost surely, we see that there is a well defined all-time minimum for the branching Bownian 
motion and we conclude that lim 2 ;^oo w(x) = 1. 

It is not hard to see by standard arguments that w must satisfy a KPP-type differential 
equation with boundary conditions: 

{q = \uj" + - uj), X > 0, 

|w(0) = 0, a;(oo) = 1. 

In fact, 0 Js{x) introduced in if finte, is solution to the same equation with the boundary 
condition t(j(0) = 0 replaced by a;s(0) = s: 


0 = gw" -I- //w' -I- I3{i0g - Us), X >0, 
Us{0) = s, Usioo) = l. 


(5) 


This is an example of the deep connection between branching Brownian motion and the KPP 
equation which goes back to McKean [16] who noticed that one can represent solutions of the 
KPP equation as expectations of functionals of branching Brownian motions. 

Until now this is very classical, however there is one unexpected difficulty here: both Q and 
(|^ admit infinitely many solutions and are not sufficient to characterize u{x). Figure gives 
several numerical solutions to Q. 



Figure 3: Solutions to Q for /3 = 1 and /x = 2. The red (bold) curve is u. 

In this work , we present three largely independent ways to characterize u{x) which are laid 
out in the three following subsections. The first approach relies on partial differential equations, 
the second gives u{x) as the expectation of a certain martingale and the third one gives u{x) as 
a power series. 

One salient property of u is that it converges to 1 rather quickly: 

Proposition 4. There exists B > 0 such that 

1 — u{x) ~ - 2 / 3 )a; large x. (6) 

Furthermore, cu(a:) is the only solution of Q which remains in [0,1) and converges that fast to 1. 
Similarly, for any s G [0, sq), there exists Bg such that 

1 — Us{x) ~ Bse~ large x (7) 

where Bi = 0, Bg > 0 when s < 1 and Bg < 0 when s > 1. 
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This fast decay was actually an unexpected result for reasons explained in Section 1.2.4 


To simplify notation, we call r the exponential decay rate of 1 — u!{x) as given in (|^: 

r ■.= ^ + IJ? -2j3. (8) 


It is the largest solution of — /rr + /3 = 0. 


1.2.1 from partial differential equations 

A first way to characterize uj{x) is to track the probability that no particle got absorbed up to 
time t. Define 

u{t,x) :=F^[K{t) = 0]. (9) 

The function u : i—>■ [0,1] is increasing in x and decreasing in t and, clearly, for each x, 

u{t,x) —>■ U}{x) as t —)• oo. Furthermore, u satisfies the KPP equation with boundary conditions 


dtu = \dxxU + + I3{u^ - u), 

u(t, 0) = 0 (Vt > 0), u(0,x) = 1 (Vx > 0), 


( 10 ) 


which, by Cauchy’s Theorem has only one solution. 

Therefore, to obtain w(x), one can in principle solve (10) and take the large time limit. This 
route leads to the following Theorem: 


Theorem 5. In regime C (ii> y/^), the function oj{x) defined by (|^ is the maximal solution 
of (|4]) such that uj{x) < 1 for all x > 0. 

More generally, ujs{x) for s < 1 is the maximal solution of © that stays below 1 and iOs(x) 
for s G (l,So] is the minimal solution of ([^ that stays above 1. 


1.2.2 Martingale representation for cj 

There is an explicit probabilistic representation of the maximum standing wave w in regime C 
(p > y/^)- Recall that A/'aii(t) is the population of all the particles in the branching Brownian 
motion with no absorption and Muveit) is the population of particles alive at time t when we 
kill at 0. We now define on the same probability space a third process based on the branching 
Brownian motion in which particles that hit 0 are stopped but not removed from the system 
(they neither move nor branch). We denote by A/'iive+abs.(^) the set of particles alive at time t 
in this model. With a slight abuse of notations we continue to write Xu{t) for the positions of 
particles when u G Mive+abs. (t)- 

Let us define the following two processes: 

^iive+abs.(t) := Zau(t):= (11) 

'^^■^live+abs. (i) 

where r is the asymptotic decay of u{x) as given in ([^. Rewriting Xu{t) = Yu{t) + ^J,t it is 
clear that {l^(t),u G A/(j,ii(t)} is simply a standard branching Brownian motion with no drift. 
Therefore Z^\\ is the usual exponential martingale with parameter r associated with the branching 
Brownian motion Y. The process .^Uve+abs. is the martingale Zaii stopped on the line t ATq (i.e. 
particles are stopped at time t or when they hit 0 for the first time). It is therefore also a 
martingale. 

Lemma 6. In regime C (/a > y/^) the martingale (Ziive+abs.(i)) ^ > 0) converges almost surely 
and in to K{oo) and therefore E’^[iL(oo)] = .^Uve+abs. (0) = e“™. 
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Following the probability tilting method pioneered by m and |8] we introduce a new prob¬ 
ability measure 

— =e^"iF(oo). 


Note that since -^live+abs. is a closed martingale we have that E®[iF(oo)|-Ft] = .^iive+abs.(^)- Thus 

dQ* _ E'live+abs. {t) 

Ft 


dP^ 


■^live+abs. (0) 


Under this tilted probability measure, the law of the process is the same as the original P* law 
except for the movement and branching rate of a distinguished particle (the spine particle ^). The 
spine moves according to a Brownian motion with drift — 2/3, branches at an accelerated 

rate of 2/3 and stops (i.e. sticks and stops reproducing) upon hitting 0. 

Theorem 7. In regime C (pL> \/^), 

1 — (jj{x) = <( 


K{oo] 


Furthermore, converges to a finite constant B > 0 when x 

I - uj{x) ^ Be-^^. 

More generally, for any s G [0, sq], one has 

I — gK(oo)' 


oo and thus, as x ^ oo. 


1 - wdx) = 


K{oo) 


and the expectation Q*(-) converges to a finite positive constant as x ^ oo. 

We will see in the proof that we can give an explicit representation of the constant B which 
appears in Proposition as the expectation of K{oo)~^ under the measure (similar to 
but with the spine particle “started at infinity”). 


1.2.3 Lj{x) as a series expansion 

The function a;(x) can be understood in terms of series expansion. Let {an}n>i be the sequence 
defined by 


ai = 1, 




13 

IjFr"^ — npir + j5 


n—1 

O-jO^n-j 

i=i 


(n 


1 


n—1 


i=i 


n>2, (12) 


(recall that r was defined in (|^, that ^r"^ — ixr + fi = 0 and that r > pL > \/2)d) and the 
function defined by the series 

^z) = Y,^nZ^- ( 13 ) 

n>l 

We have 


Proposition 8. The radius of convergence IZ of ^ is non-zero and there exists B G (0,7^) such 
that 

u:{x) = 1 - $(Se-™). 

More generally, for any 0 < s < sq, there exists a number Bg such that 

0 Js{x) = 1 — ^(Bse~^^) for all x > 0 such that \Bs\e~'^^ < IZ. (14) 

s ^ Bg is decreasing, positive for s < 1, zero for s = 1, and negative for s > 1. In particular, 
for s < 1, the condition \Bg\e~'^^ < IZ is automatically fulfilled. 
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The proof is contained in Section 2.2.4 Numerically, it seems that TZ is large enough that 
\Bs\e~^^ < TZ for all s G [0, sq] and all x > 0, but we haven’t proved that point. The represen¬ 
tation (14) makes it very easy to compute numerically iOg by first computing d>( 2 ;), see Figure]^ 
the value sq is then obtained as 1 minus the first minimum of for negative arguments. This 
follows easily from the facts that w' = 0, w” < 0 and that rc^(O) < 0 for all s G (1, sq). 



Figure 4: The function z i—>■ d>( 2 ;) for /3 = 1 and /x = \/2 (left) or /x = 3 (right). Numerically, one 
has So ~ 1.35, « -0.859, Bq ^ 0.564, « 3.14 for /x = \/2 and sq ~ 14.11, Bs^ « -39.86, 

Bo « 0.969, n PS 72.8 for /x = 3. 


1.2.4 About the asymptotic decay (|^ 

The KPP partial differential equation, 

dth=]^dxxh +(3{h^-h), (15) 

where 0 < h{t,x) < 1, h{t, —oo) = 0 and h{t, -|-oo) = 1, describes how a stable phase (/i = 0 on 
the left) invades an unstable phase (h = 1 on the right). It is well known that it admits travelling 
wave solutions of the form 


h{t,x) = hfj^{x - fit), 0 </x^ < 1, hfj_{+oo) = 1, /x^(-cx)) = 0, 

for any velocity /x greater or equal to \/2)0. The travelling wave x i—)• h^{x) is then solution to 


+ V) = 0, h^{-oo) = 0, h^{+oo) = l. 


(16) 


The solution to (16) is unique up to translation. 

Equation (j^) for is very similar to equation (|^ for uj{x), but h^, surprisingly, does not 
have the same asymptotic behaviour as uj for large x, in the region where and w are close 
to 1. Indeed, linearising (16) around 1, one gets 


^(1 - hfj,)" + /x(l - h/j,)' -h /3(1 - hii) = 0, [linearized] 


(17) 


(a term of order (1 — h^)^ has been neglected) and the general solution to 0 is, for some 
constants A and B, 


1 - hfj_{x) 


(Ax + B)e~'^^ 
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for/X > \/2^, ... . . 

V / ; [linearized] (18) 

for /X = y/^. 




















For X large, hn is close to hn with the meaning that for some constant A and B, 1 — /i„ ~ 1 — h. 


(see the discussion in Section 2.2.3). Of course, if A ^ 0, the term in factor of B is negligible 
compared to the term in factor of A and the A term alone is an equivalent to 1 — h^. When 


solving (16), it turns out that the solution has a non-zero A term and that, therefore. 


h^{x) 


for 

\Axe~'^'^ for fi = -v/2^. 


(19) 


where A depends on /x. 

We now consider equation (Q for uj(x). Of course, the boundary condition of (Q is not 
sufficient to determine a unique solution, and for a range of values of c there exists a solution to 


0 = -v" + fj,v' + f3{v^ — v), u(-|-oo) = 1, u(0) = 0, v'{0) = c. 


( 20 ) 


(The difference with equation (Q for u}{x) is the added condition u'(0) = c. Figure shows 
several solutions.) 

One can then do, as above, a large x analysis of v and, the partial differential equation being 
the same, one finds again that 1 — u ~ 1 — {x large) as given in (18) for some c and ^ 


dependent values of A and B. Generically, A is non-zero and 1 — v decays as 1 — in ( [19[ ) (up 
to a multiplicative constant; the A is usually different and can even be negative, see Figure]^. 

However, for a well chosen value of c (depending on /x), one has A = 0, and the asymptotic 
decay of \—v is given by the B term (that is: it decays much faster). The meaning of Proposition]^ 
is that uj{x) is precisely that very special solution to (|^ that decays unlike all the other ones 
and unlike the travelling wave h^. 

It is also interesting to remark that an equation very similar to (Q appears in the study 
UDE] of the the extinction probability of a branching Brownian motion with absorption and 
supercritical drift /x > —y^2)0 (regimes B and C): let 0{x) = P*[A/'iive(c>o) = 0] be the extinction 
probability when the system is started from x; then one has 


'0 = f0" + ^0' + /3(02-0), x>0, 

0 (0) = 1, 0(oo) = 0. 


( 21 ) 


Equations (Q and (21) differ only by their boundary conditions; however (21) has a unique 
solution, whereas (|^ has many. 

A possible way to understand the difference is that an asymptotic analysis of 6{x) for large x 
similar to (18) yields only one possible exponential decay: 9{x) Aexp[—(;U -k /x^ -k 2/3)x] for 


some constant A, which means that, up to translations, there is only one solution which does 
converge to zero at infinity, whereas for uj{x) there were two possible exponential decays and 
infinitely many solutions. Otherwise said, if one were to impose 0(0) = 1 and 0^(0) = —c, there 
would only be one value of c for which 0 would converge to zero at infinity. 


1.3 A travelling wave result 

In regimes A and B (^ < \/2)d) one has oj{x) = 0 because K{oo) > 0 almost surely, which is not 
very interesting. What is more interesting is the way that u{t, x), defined in ([^ as F^[K{t) = 0], 
converges to zero: it does so by assuming the shape of the critical travelling wave of the KPP 
equation. Let us recall quickly the well known facts on this critical travelling wave. 


Consider the KPP equation (15) without drift on the whole line with Heaviside initial con¬ 
ditions: 

f dth = \dxxh + I3{h‘^ - /i), 

[ h(0, x) = 0 (Vx < 0), h(0, x) = 1 (Vx > 0). 


( 22 ) 
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It is well known that h{t, x) is the probability that the leftmost particle at time t of a branching 
Brownian motion started at x is to the right of zero. Furthermore, this probability converges to 
the critical travelling wave in the following sense: 


h{t, rrit + x) ^ /i*(x) uniformly in x as t —)• oo 


with 


m : 


■■= Vwt- 


2 ^ 2 ^ 


log t + Cste 


(23) 


and where := h is the travelling wave moving at the minimal possible velocity \/ 2 ^, 
see (16). To fix the invariance by translation, we impose the further condition /i*(0) = 1/2: 


'Q = \K + ^K + p{hl-K), 

/i*(—oo) = 0 , /i*( 0 ) = 5 , /i*(+cx)) = 1 , 


(24) 


and the solution to (24) is now unique. 


Adding a drift +^dxh to (22) would only shift the solutions by —fxt and would make (22) 


very similar to ( | 10 [ ): the only difference would be that u is defined on M'*' and h on M, but as 
both equations converge quickly to zero around the origin in regimes A and B (^ < \/^), this 
difference turns out to be minimal and one has 

Theorem 9. In regimes A and B (fi < y/^), there exists a constant C depending on p, and /3 
such that 

u{t, X + mt — pt + C) ^ h^: (x) uniformly in x as t —)• 00 


where mt is given by (23) and /i* is the solution to (24). 


It is interesting to compare this result about the behaviour of u{t,x) = ¥^(K{t) = 0) when 
p < \/2)5 to the behaviour of the extinction probability u{t,x) = P^(A/'iive(i) = 0) when p < 


— \/2^. It is not hard to see that u satisfies the same equation (10) as u with different boundary 


conditions, which is 1 minus the boundary condition in ( | 10 [ ); namely u solves 

j dtu = \dxxU + pdxU + - ft), 

[ u{t, 0) = 1 (Vt > 0), tt(0, x) = 0 (Vx > 0), 


(25) 


What is particularly striking is that in the critical case p = it is known since Kesten 

m (see also 0) that to survive up to time t one must start with an initial particle at position 
X = This means that if u{t,x + mt) converges to some limit front shape then the centering 

term giving the position of the front fht has to be of order ct^^^. However the convergence of the 


solution of (25) to a travelling wave is at present an open problem. 


2 Proofs 


The proofs section are presented mostly in the same order as the results: Section 2.1 contains 
the proofs to Theorems an d about the structure of the functions ujs{x) = 


m 


regime C {p > \/2^). Section 2.2 contains the proofs to Proposition]^ Theorem]^ Lemma 


Theorem]^ and Proposition]^ about the different representations of a;(x) = wo(x) = P*[iL(oo) = 
0] in regime C {p > \/2)5). Finally, Section 2.3 contains the proof of Theorem ]^ about the 
establishment of a travelling for u{t,x) = ¥^[K{t) = 0] in regimes A and B (// < -v/2/3) and the 
asymptotic behavior of sq established in Proposition ]^ is proven last. 
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2.1 Proof concerning the tail behaviour of K{oo) 

In this section we consider exclusively regime C (/i > \/2^) and we focus on the problem of 
the exponential moments of K{oo). We first establish some properties of uJs(x) = as 

defined in Q and proceed to prove Theorems]^ and Proposition]^ We then prove the asymptotic 
behaviour (??) to complete the proof of Theorem]^ 

The first property we need is that for a given s, the quantity 0 Js{x) is either finite for all x > 0 
or infinite for all x > 0 : 


Lemma 10. For a given s, 

(3x > 0 : u}s{x) < +oo) (Vx >0 : uis{x) < +oo). 

Proof. Fix s > 0, X > 0 and y > 0. There is a positive probability, which we note €{x,y), that 
the initial particle starting from x reaches position y before any branching or killing happens. 
Then 

a;,(x) = > e(x,= e{x,y)uM- 

Therefore, if a;s(x) is finite, then 0Js{y) is also finite. □ 

Remark: in the following, we write < oo when the conditions of the lemma are met. 
Clearly, this is the case when s < 1. Furthermore, as s i—)> 0 Js{x) is obviously increasing, if 
oJsQ < oo for some sQ) then ovg < oo for all s < sq. 

When cjs < oo, it is clear by standard arguments that a;s(x) is solution to 


0 = 5 ^'' + //W' + - LVs), 

UJs{0) = s. 


(26) 


Let us now prove Theorem]^ A slightly more general result is given by the following Lemma: 
Lemma 11. If eog < oo then, for any x > 0 and h>0, 


UJs{x + h) = U}uJs{h){x)- 

Setting s = 0 and renaming u}o{h) as s gives the first line of the Theorem. Once we have 
proved that sq exists, setting s = sq and renaming ujsg(h) as s gives the second line of the 
Theorem. 


Proof. Instead of starting our branching process at position x and killing particles at 0, it is 
here more convenient to think of the process as started at 0 and particles being absorbed at —x. 
This allows to couple different values of the killing position. In particular, if Tlx designates the 
particles stopped when they first hit —x and Kx{oo) is the number of particles in Hx, we have 
that 

ujs{x + h) = E[ =IE[ ujs{h)] =w^4fe)(x), 

where Kj^\oo) is the total number of descendent of the particle u which are killed at —x — h 
(which by translation invariance of the branching Brownian motion and the branching property 
is an independent copy of K^ioo)). □ 

We now have a monotonicity result: 

Lemma 12. 

1. // s < 1; X I—)• 0Js{x) is an increasing function converging to 1. 

2. // s > 1 and uig < oo, x i—>■ a;s(x) is a decreasing function converging to 1. 
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Proof. Once the incr easing/deer easing part is proved, the fact that the limit is 1 is obvious: from 
its definition, it is clear that a;<j < 1 if s < 1 and > 1 if s > 1. Assuming uig is increasing or 
decreasing (depending on s), it must have a limit, and from (26) that limit must be 1. 

From its interpretation as the distribution of the all-time minimum of a branching Brownian 
motion, see (|^, it is furthermore clear that wq = w is an increasing function. Then, the coupling 
provided by Lemma 11 (or more simply Theorem]^ implies that iOg is an increasing function for 
all s < 1. 

Therefore, it only remains to prove that for s > 1, is decreasing when it is finite. 

Assume s > 1 and < oo. We first show that Ug is monotonous by considering two cases: 


If ca((0) > 0 then, for all /i > 0 small enough, u!s{h) > s. But, for x fixed, s i—?• lOs(x) is a 
strictly increasing function so co^^(h)(x) > ^^(x). Then by Lemma 
for all X and all /i > 0 small enough: Wg is increasing. 
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LUg(x + h) > Ldg{x) 


• If w(,(0) < 0 then, for all /i > 0 small enough, ojg{h) < s because in the limit case w((0) = 0, 
one has w/(0) < 0 from (26). Then, as in previous case, ujg{x + h) = a;^^(/i)(x) < ujg{x) for 
all X and all /i > 0 small enough: ojg is decreasing. 

It now remains to rule out the possibility that ojg is increasing for s > 1. Imagine that s > 1 
and (jjg increases. Then, from (26), oj”{x) < —2/3(Wg(x) — a;s(x)) < —2/3(s^ — s) and w((x) < 
a;((0) — 2/3(s^ — s)x, which becomes negative for x large enough, in contradiction with the fact 
that ujg increases. So ujg must decrease for s > 1. □ 

We need now to characterize the values of s for which < oo. 

Lemma 13. Assume s > 1. If there exists a function v which solves 

-v” 

2 ^ 

v{0) = s, 


0 = Iv" + ^iv' + /3{v‘^ - v), 


v{x) > 1 (Vx > 0), 


(27) 


then ujg < oo. 


Remark: The converse is obvious: when ojg is finite, it is one of the solutions to (27). This 
Lemma allows to define 


So = sup {s > 1 : cus < oo}, 

= sup {s > 1 : a solution to (27) exists}, 

and because s i— ?• ojg{x) increases, one has ojg < oo for all s < sq- 

Proof. We present two proofs: one probabilistic and one analytical. 


(28) 


Choose s > 1 such that (27) has a solution v. We introduce the process 


Mt := H 

^£.^live+abs. (L 


where we recall that A/iive+abs. (i) is the set of particles in the branching Brownian motion where 
particles are frozen at the origin, see Section |1.2.2[ 

Mt is a positive local martingale and therefore a positive super-martingale which thus con¬ 
verges almost surely to Moo- Observe that under 


t^(a:) = Mo > 


But since for all t > 0 one has 


Mt > x(0)^W = 
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we see that M^o > and therefore 

w,(x) <r;(x) <oo. 

The same result can be proved analytically through the maximum principle. Let us introduce 

Usit,x):=E^[s^^% (29) 


which is clearly solution to 


dtUs = ^dxxUs + fJ-dxUs +/3{u^ - Us), x>0, 
Us{t,0) = s {'it> 0), tis(0,x) = 1 (Vx > 0). 


(30) 


(Compare to .) With V as above, one clearly has Vx > 0, ris(0,x) < v{x). Therefore, by the 
maximum principle we have that 


Us{t, x) < v{x), Vt > 0, X > 0 


and as Us{t, x) uJsix) as f —>■ oo we see that ujs{x) < v{x) < oo. □ 

It is obvious that sq defined in ( |28[ ) depends only on the ratio by a simple scaling 

argument: the branching Brownian motion with drift /r and branching rate 13 is transformed, 
when time is scaled by A and space by ^/X, into a branching Brownian motion with drift /U\/A 
and branching rate (3\. In particular, ujs,p,ij,{x) = i3\ ~ V~P) with the 

obvious new notation. What remains to be shown are the following properties of LOg' sq is finite, 
iOgQ is finite, a;(|j(0) = 0 and so > 1. 

Lemma 14. sq < oo, that is: K{oo) does not have exponential moments of all orders. 

Proof. For the system started from x > 0, consider the following family of events for n G N: 


.4 


n 


K{oo) = n, and 

< for all integers i < n, K{i) = i, and 

for all integers i < n, there is at time i only one particle alive and it sits in [x, x + 1]. 


In words, for each i G {0,1,..., n — 1} there is one particle alive at time i, it sits in [x, x + 1], 
and during a time interval one, this particle splits exactly once, one the offspring gets absorbed 
and the other is again in [x,x + 1] at time i + 1. The one particle alive at time n generates a 
tree drifting to infinity with no more absorbed particles. 

Let ey^zdiZ be the probability that a particle sitting at y has, during a time interval one, 
exactly one splitting event with one offspring being absorbed and the other one ending up in dz. 
Define furthermore 

rx+l 

q = min / ^y,zdiz. 

y&[x,x+l\ Jx 

q is the minimal probability for a particle sitting somewhere in [x, x + 1] to have, during a time 
interval one, exactly one splitting event with one offspring being absorbed and the other one 
ending up in [x,x + 1]. It is clear that q > 0 and that, furthermore, 

r^{K{oo) = n)> F^{An) > g”P^(iL(oo) = 0). 

This implies that = oo and that sq < l/<? < oo. □ 

Lemma 15. uig^ < oo. 
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Proof. If So = 1, this is trivial as oji = 1. Assume now sq > 1 and let us fix x > 0. For any 
1 < s < So, as is decreasing, one has a;s(a;) < s < sq- This implies that (using the monotone 
convergence Theorem) Wsq(x) = lim^^sp ujs{x) is finite which entails the result. □ 

Lemma 16. = 0. 

Proof. We already know that ^(^^(0) < 0. Assuming w(q( 0) < 0, one could continue the function 
ujso to negative arguments using ( |27| ) and one could find a Xo < 0 such that ^^^(xo) > sq; then 
the function x i—)> a;sp(xo + x) satisfies ( |27[ ) with s > sq, which is a contradiction. □ 

The proof that so > 1 for all ^ > \/2P is divided into two steps. First, we show that sq > 1 
in the critical case /U = \/2)d. Then we conclude by proving Proposition!^ which states that sq 
is an increasing function of /x/y(d. 


Lemma 17. In the critical case ^ for s > 1 small enough, there exists solutions to 

that is So > 1. 


(271, 


Proof. Assume /r = \/2)d. After the change of variables i{x) := e^^{v{x) — l), (|27[) reads 


1 

2 


i" + = 0 . 


(31) 


Let us consider the solution to (31) with £(0) = T(0) = e for some e > 0. We want to 
that Vx,£(x) > 0 if e is small enough. Assume otherwise and call xq = inf{x > 0 : i{x) 
Then, as i"{x) < 0, we have i{x) < e + xe and thus on [0,xo] (where f'(x) > 0), 


prove 
= 0 }. 


/'(x) > -2l3e^{l + x)‘^e-^^^. 


We conclude that 


t-xo poo 

£'{xo) > e - 2pe‘^ / (1 + xfe->^^ dx > e - 2pe^ / (1 + xfe->^^ dx, 

JO JO 


which is strictly positive for e small enough. This contradicts the definition of xq and thus we 
have found a solution of (31) such that £{x) > 0 for all x > 0. Then v{x) = 1 + £{x)e~^^ is a 
solution to (27) started from s = 1 + e; in other words sq > 1 + e in the /r = \/2(5 case. □ 


We now proceed to prove that sq is a strictly increasing function of ^/y/p. 


Proof of Proposition Let us fix /i > y/2^ > y/2j32. One can easily construct two branching 
Brownian motions with parameters {p., Pi) and {p,, P 2 ) on the same probability space to realise 
a coupling so that the particles of the second one are a subset of the particles of the first one. It 
is then clear that for any s > 1 one has 




(32) 


(with the obvious extension of notation) so that 

soA/AA) < soiix/y/Ai)- 


(33) 


This already gives non-strict monotonicity and concludes the proof that sq > 1 for all p., /? with 

h > VW- 


We can now prove that the inequality (33) is strict. Assume otherwise; one would have 

= tJso,/ 32 ,/i( 0 ) = 'So (where sq > 1 would be the common value), ' ■' 

2 


u 


UJ 


j(0) = 0 (from Lemma 16) and, from (26), ^ 




(0) = -/3i(sg - so) < 


( 0 ) = 


so,132,U _ 

—P 2 {sq — So), which would imply by Taylor expansion that for x > 0 small enough a;sp_^^^^(x) < 

oJso,I 32 ,uA) contradiction with (32). □ 
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Finally, the only remaining point to complete the proof of Theorem is the asymptotic 
behaviour Q, i.e. the assertion that for x > 0 fixed 


Pn{x) := P^[ii:(oo) = n] 


-Wsn(a;) 


28^72 2 Y^7r/3(so - 1) 


as n —)■ oo. 


(34) 


Write D{z, r) for the open disc of the complex plane with center 2 ; G C and radius r. We 
extend the definition of s 1 —?• 0 Js{x) to s G C: 


a;s(3:) = ^ s^pn{x), Ms G C. 


(35) 


n>0 


This quantity is analytical on D{0, sq) because the Pn in (35) are positive and the first singularity 
on the real axis is at sq- Furthermore, it is finite on -D(0, sq) by uniform convergence because it 
is finite at sq- 

The arguments we use are extremely close to those used by Maillard in |15) . The key ar¬ 
gument, which improves on usual Tauberian theorems, is an application of |9l Corollary VI.1] 
relying on the analysis of generating functions near their singular points. We need to show that 

Lemma 18. Fix x > 0. There exists > 0 such that s 1 —>■ uJs(x) is analytical inV = D{sq, Vx) \ 
[so, 00 ), and 


dscas(x) 


2\//3(so - s){sl - So) 


as s ^ So, s & V, 


(36) 


and that 


Lemma 19. Fix x > 0. There exists e > 0 such that s 1 —)■ a;s(x) is analytical on 11(0, sq + e) \ 
[so, 00 ). 


Then, applying Corollary VI. 1 in [9] to s 1 —)• dsLOs{x), Lemmas 18 and 19 lead to 


(n l)pn+i{x) 


2so+^n2 y^7r/3(so - 1)’ 


as n —>■ 00 , 


(37) 


which obviously implies (34). 


Proof of Lemma We know that a;(jj(0) = 0 and = —2/3(sq — sq) < 0. Since Wsq solves 

the KPP traveling wave differential equation, for each x > 0 we can extend z 1 —>• ^^^(x -|- z) 
analytically on a neighborhood of zero in C (see e.g. |18)1. In particular for x = 0 we have the 
following expansion: 


Usoiz) = So + 


^so(6)^2 I „/'^2 


z^ + o{z^) as z —0. 


(38) 


The function is analytic and zero is a zero of order two of ujsg(z) — sq, by Theorem 10.32 of 
|18| there exists ri > 0 and a function ip analytic and invertible on D{0, ri) such that 


UJsoiz) = So + 




This means that 


z = Ip 


-1 ^soiz) - So] _ So-Usgiz) 


<(0)/2 


= Ip- 


Ksl - So) 


(39) 


(40) 


for any z in H(0,ri) such that lOso{z) 0 (so, 00 ) (so that the right-hand side is well defined and 
is analytic on this domain when using the standard definition of the complex square root). 


Recall from Lemma 11 that for any non-negative real x and z one has 

^^LOsoiz)ix) =UJso{z + x). 


(41) 
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Replace the z in the right-hand side by its expression (40) and write cOsoi^) s to obtain 


a;s(x) = (XJ so [i) 


-1 


/ So - g 
^(Sq - ■®o) 


-|- X 


(42) 


for s G a;sQ([0,ri)) = (sq — 5 " 2 )'So] for some r 2 > 0. But (42) is an equality between analytical 
functions as long as s G D{so,rx) \ [so,oo) for some > 0 small enough (one must have 
D{so,rx) C WsQ (i4(0, ri)) for to be analytical, which is possible by the open mapping 

Theorem, and one must have ■ ■) small enough for Usq to be also analytical). From the 

analytical continuation principle, (42) must hold on the whole D(so,rx) \ [sqjOo) domain. Now 
differentiate with respect to s to get 


(-0 


-U/ 


dsOJsix) = - 


SQ-S 

/3(so-so) 


2V/3(so - 'S)(so - ■So) 


U) 


so 


■0 


-1 


I Sq- S 

P(4 - ■So) 


-|- X 


(43) 


yielding 


dgUJsix) - -- 


(r')'(o) 


uj'(x) as s So in B(so,rx) \ [so,oo). 


2\//3(so - s)(s^ - So) 

A straightforward computation shows that (V’~^)^(0) = 1, which concludes the proof. 


(44) 


□ 

Proof of Lemma [7^ 0 Js{x) is already analytical on T>(0, sq)- To prove the Lemma it is sufficient 
to show that it can be analytically extended around any point s G dD(0, sq) \ {so}- Indeed, by 
the finite covering property of compacts one can then show analycity on a open containing the 


compact dD{0,so) \ D{so,rx/2) with defined in Lemma 18 and then we conclude with the 
help of Lemma 

So it now remains to see why s i— ujs{x) can be analytically extended to neighboorhoods of 
any s 7 ^ sq with |s| = sq. This is essentially the content of Lemma 6.2 in |15) . As in Maillard, 
we define 

a{s) := w((0), (45) 

where we recall that the prime is a derivatice with respect to x. We first show analycity of a(s) 
on D{0, So) by writing an integral representation of a(s): multiply ([^ by exp [{fj, — _|_ 2/3)x] 

and integrate on x G [0, 00 ). Integrate several times by part to get rid of the derivatives of LOg] 
one is left with 


i{s) = (// - V^2 + 2/3)s + 2/3 r dxujs{xfei>^-'^^^)^. 


(46) 


For any x > 0 and s G D{0, sq) one has |ws(x)| < cOsqIx) < sg- This implies that the convergence 
for X close to infinity of the integral in ( |46[ ) is uniform on the disk s G B(0, sq}- As s i — cas(x) 
is analytical on B(0, sg), this is sufficient to ensure that s 1 —?• a(s) is also analytical on B(0, so). 


Furthermore, notice that the series (35) defining uis{x) converges uniformly on s G D{0,so) 
because we know it converges absolutely (all the Pn{x) are non-negative) at s = So. This implies 


that s I—)• a;<j(x) is continuous on L)(0, sq ) and, from the expression (46), so is s 1 —)• o(s) (by 
dominated convergence Theorem since |a;s(a;)| < so on the closed disc). 

We proceed to show that a(s) can be extended analytically around any point s 7 ^ So with 
|s| = So and show that the property extends to 0 Js{x). 


In Lemma 11 we showed for any s G [0, sq ] and any x > 0 and h > 0 one had 


uJs{x + h)= uj^^(^h){x). 


(47) 
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One can check that the proof of Lemma 11 extends to complex s so that (47) remains valid for 
any s £ C such that u}s{x) is hnite. 


For hxed (complex) s, by deriving ( |47[ ) with respect to h and then setting h = 0, one gets 

Ug{x) = a{s)dsUs{x). (48) 

Derive again with respect to x, and then set x = 0: 

a;"(0) = a{s)dsa{s), (49) 


so that the differential equation (26) on x i—)• u}s{x) applied at x = 0 gives 

0 = ^a(s)9sa(s) + /xa(s) + /5(s^ - s), 


(50) 


which is equation (3.5) in |15| . This equation is valid for all s £ D{0, sq)- 

We now use the Fact 5.2 in m, which we reproduce here with some change of notation for 
clarity: 

Fact 5.2 in m Let 77 be a region in C and s i—>■ <j){s) analytic in H. Let G be a 
region in such that {4>{s), s) £ G for each s £ H and suppose that there exists an 
analytic function / : G —)• C such that 

(p'{s) = f{(j){s),s), Vs G 77. 

Let s* G 977. Suppose is continuous at s* and that ((p{s*),s*) £ G. Then s* is 
a regular point of 4>{s) i.e. (p{s) admits an analytic continuation at s*. 


We apply to our case with (f) = a, H = 77(0, sq) \ {1} and G = C* x C. From ( [^ , the only 
candidate values of s in 77(0, sq) such that a(s) = 0 are 0 and 1, and we know that a(0) > 0, so the 
condition “(a(s),s) G G for each s G 77” is verihed. The function /(a, s) is obtained from ( |50[ ): 
/(a, s) = — 2 ^ + 2(s — s^)/a, and is obviously analytical on G; we have already shown that a(s) 
is continuous on 77(0, sq). Therefore, for any point s* £ 977(0, sq) such that o(s*) 7 ^ 0 (because 
we want (a(s*), s*) G G), the function a(s) admits an analytic continuation at s*. We know that 
a(so) = 0 , and we prove now that one has a(s*) 7^ 0 for any s* £ 977(0, sq) \ {so}, which will 
conclude the proof that o(s) can be analytically continued around any point in 977(0, sq) \ {sq}- 
From (35) one can write a{s) as a series: 


a{s) = ^p;(0)s”. 


(51) 


We know that Pi(0) < 0 (because pi(0) = 1 and pi(x > 0) < 1) and that for n 7 ^ 1, p'„(0) > 0 
(because Pn(0) = 0 and Pn{x > 0) > 0). Since a(so) = 0, we write 


rif^l 


(52) 


All the terms on the left hand side are non-negative and inhnitely many of them are non-zero 
since (50) does not have polynomial solutions. Thus, for any s* £ 977(0, sq) \ {sq} one has 


n^l 


< 5^K(o)«o- 

riT^l 


(53) 


In particular, a{s*) 7 ^ 0 because it is the sum of two terms (the term n = 1) with 

different moduli and s 1 —)• a(s) can be extended analytically around s*. 
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We now show how the analycity of a(s) translates into analycity of a;s(x). First derive (47) 
again but this time with respect to x, then set x = 0, and rename h into x to obtain 


ujg{x) = a[a;s(x)] = a{s)dsU}s(x) 


(54) 


where we used (48) for the second equality. 


For each given s* G dD{0, so)\{so} we consider a neibourhood V of s* where a{s) is analytical 
and we apply again Fact 5.2 in |15| to prove that s i—?• a;s(x) is also analytical around s*. This 
time, we take </>(s) = iOs{x) and /(w, s) = a{uj)/a{s) from ( [54| ). We pick H = D{0, sq) \ {1} and 
G = D{0,so) X {H U V). For any s G 74(0, sq) \ {so} one has |a;s(x)| < Ws(,(x) < sq so that the 
condition “((/>(s),s) G G for each s G 77” is satisfied. We have already shown that s i— ijJs{x) is 
continuous on 74(0, sq), so we conclude that s* is a regular point of s i—)• a;s(a:). 

□ 


2.2 Proofs concerning uj{x) 

in this section we consider exclusively regime C (/r > -v/2^) and prove our results concerning the 
properties of u and Ug- 


2.2.1 Proof of Theorem [5] 

We need to prove that uj is the maximal solution remaining below 1 of the differential equation (Q . 
This is an elementary application of the maximum principle again. Suppose that v is any solution 
of (Q which stays below 1. Since r; is a standing wave solution of (10), that is u{t, x) = v{x) for 
all t > 0 is a solution of 

dtu = -dxxU + gLdxU + / 3 ( m ^ - u), 


and since v(x) < 1 for all x > 0 we have that 

v(x) < u(t, x), 


Vt > 0,Vx > 0 


where u is the solution of (10). As for each x we know that 1 1—)• u{t, x) uj{x) we conclude that 
uj{x) > v{x) and therefore ui is the maximal solution of (Q bounded by 1. The same argument 
is easily generalized to the case of an arbitrary value of s G [0, sq]. 


2.2.2 Proof of the martingale representation, Theorem 

We start by proving Lemma i.e. that the martingale (.^iive+abs.(^)T ^ 0) converges P-almost 
surely and in to K{oo) and therefore that E^[7L(oo)] = 

Proof of Lemma 0 Recall that by ( [TT| ) 

^live+abs.(i) := ^ ^-rX„{t)^ Zall(t) := ^ ^-rXu{t) ^ 1-55^ 

^e.^live+abs. (^) 

are positive martingales which therefore converge P-almost surely to their respective limits 
Ziive+abs. -^all- Furthermore, as Z^\\[f) is the usual additive martingale with parameter 
r>V^, one has Z^n = 0. As the bounds 

7L(t) < Ziive+abs.(i) < 7L(t) -|- Zaii(t) 

always hold, it is clear that Znve+abs. = K{oo). The only thing left is to show that the convergence 
also holds in L^. 
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We start by recalling the description of the measure which is defined by 

^ ^live+abs .(0 

dP^ Tt ~ Zu,e+abs.( 0 )' 

Standard arguments (see [10]) allow us to conclude that under the process behaves as follows: 
for t > 0, there is a distinguished line of descent (the spine) denoted ^{t) G A/iive+abs. (^) ■ Under 
the particle ^ moves according to a Brownian motion with drift — 2/3 and therefore 

almost surely hits 0 in finite time; we call rg = inf{t > 0 : = 0} the time at which it 

reaches 0. For t < r^, the spine branches at rate 2/3 creating non-spine particles which start 
new independent branching Brownian motion behaving according to the usual P law. After 
the spine particle is frozen at zero (no motion, no branching). Observe that is actually the 
projection of the measure just described since under we do not know which is the spine 
particle 

To prove the convergence of Ziive+abs. (0 towards its limit Zuve+abs. := limt ■^live+abs.(i), it 
is sufficient to show that 

Q (-^live+abs. < Oo) = 1 . 

As the time at which the spine is absorbed at 0 is Q-almost surely finite, there are only 
finitely many branching events from the spine Q-almost surely as well. At each of these events, a 
non-spine particle u starts its own independent P branching Brownian motion and we call Ku{oo) 
the total number of particles frozen at 0 that are descended from u. Let us also call ■^ijve_|_abs 
the analogue of the limit Ziive+abs. (but we sum only on particles descended from u) and is 
the same as -^b^e+abs without any absorption or freezing at 0. It is clear as above that 

K.(oo) < 4;Lb.. < + zS'- 

and that = 0. We conclude that ~ ^u{oo) < oo, Q-almost surely, and finally 

Ziive+abs. = K{oo) < OO Q-almost surely. 

Observe that since K(oo) > 1 almost surely under Q, we have Q ~ P. Thus we know that under 
P Ziive+abs.(i) -t Ziive+abs. = K{oo) in Lb Hence, E^(L:(oo)) = Ziive+abs.(0) = □ 


We now move to the proof of Theorem 


Lemma 20. Recall ui{x) 
Then 


P^(iL(oo) = 0) and lJs(x) := for 0 < s < sq as usual. 


1 - io{x) = 


- 1 
-K{oo) 


e 


—rx 


and 


1 - Ws(a;) 




ri _ s^(oo)-| 

^ ^ p-rx 

- K{oo) \ 


Proof. As K{oo) > 0 Q*-almost surely, it is sufficient to prove the second assertion. Using that 
K{oo) = .^live-i-abs. P*-almost surely. 


1 - = P-ll - = P-Kl - 


= P^[(1 - 

= Q"[ 


•^live-l-abs. ( 0 ) ■^live+abs. 


X f ^live+abs. ( 0 ) 




live+abs. 


■^live+abs. -^live+abs. ( 0 ) 

(1 - 


{-^live+abs. ^0 


}] 


= Q"[ 


I _ gK{oo) _ 

K{oo) - 


□ 
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Since we already know from Proposition (its proof is analytical and is independent of the 
present discussion) that (1 — U 0 s{x))e'^^ tends to a constant -B > 0 , it is now clear that the Q* 
expectations in Lemma also converge to B as x —oo. However, we are now going to define 
as the law of the process under which we can couple all the Q* together and interpret the 
limit constant B as the expectation of a limit variable under Loosely speaking, we want 

to be the law of the process where the spine particle starts at x = +oo before drifting to 0 . 
In fact it is easier to reverse time and have the spine start at 0 and drift to +oo. 

We start with the following Lemma: 

Lemma 21. Let (Y{t),t > 0) be a Brownian motion with drift +y/ — 2/3, started from 0 
conditioned to never hit 0. Otherwise said, Y is solution of the following stochastic differential 
equation 

j dy (t) = dBt + _ 2/3 coth — 2/3 H (t))dt if fa > 

|dy (t) = dBt + if h = VW- 

Let (ti) be a Poisson point process on with intensity 2/3. For each i > 0 start a branching 
Brownian motion with law and call Kt the total number of absorbed particles at 0 for this 

process. 

Fix X > 0, then the distribution of the variable K{oo) under is the same as that of 

K^ioo) := 1 + Yj 

under where Tx ■= sup^>o{y(t) = x}. 

This result should be clear once it is realized that the process Y is the reversed path of the 
spine 

Proof. We only treat the case fi > \/^ since the zero-drift case is similar. The only thing we 
need to prove here is that if < r^) is a Brownian motion with drift — /i^ — 2/3 started 

form X and stopped at time := inf{t : f^(t) = 0}, then 

{ifit),t < = {{Y^Tx -t),t< Tx),Tx}. 

This follows for instance from |^. 

□ 


The upshot of Lemma 21 is that we can now construct the variables K(oo) under for 
all values of x simultaneously. We write Q°° for the joint law of the variables {(Y(t),t > 
0), (L)igN, Ki) described above. Then under clearly (iL^(oo), x > 0) is an increasing process 
in X. We call K°°{oo) its limit which is also the total number of particles absorbed at 0 under 


Lemma 22. We have that K°°{oo) < oo -almost surely. 

Proof. We start with the fa > \/2^ case. First observe that for e > 0 fixed, there exists almost 
surely a random zq G N such that 

\/i>io, Y{ti) > (—-e)»- 

This simply comes form the fact that Y{t)/t —>• — 2/3 and ti/i —)■ (2/3)“^ almost surely. 

Now, 1 — oj{x) < because under we have that Ziive+abs. > 1 almost surely and 

therefore Q^(l/yiive+abs.) < 1- Hence, for any z > zq we have that 

Q^{K^ > 0 ) < 1 - 
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for some positive constant c. Thus a straightforward application of Borel-Cantelli Lemma shows 
that almost surely, there exists jo G N such that Vi > jo, iLj = 0, which yields the desired result. 
The zero-drift case is similar. One just needs to start the argument by observing that for e > 0 
fixed, there exists almost surely a random zq G N such that 

Vi > io, Y(ti) > 


where c is a constant. The proof then follows as before. 

The following Lemma completes the proof of Theorem 
Lemma 23. For s G [0, sq] small enough, we have that 


□ 


gK{oo) 

K{oo) 


gK°°(oo) 

K°°{oo) J 


(56) 


Proof. The monotone convergence Theorem applies when s < 1 so we suppose 1 < s < sq- 
Observe that the map 1 1—>■ s^/t is decreasing on [1,1/logs] and increasing on [l/logs,oo). Thus 
we write 


gK{oc) 

K{oo) 


gK-{oo) 

K^{oo) 

gK-ioo) 

K^{oo) 

„X°°(oo) 


K=={oo) 




K°°{oo) 




K°°(oo) 

gK°°(oo) 

K°^{oo) 


where the first convergence comes from the dominated convergence Theorem and the second 
from the monotone convergence Theorem. 

□ 


2.2.3 Proof of Proposition 

We consider here all the solutions x i— v{x) to Q that remains in [0,1). By Cauchy’s theorem, 
a solution to Q is entirely determined once the derivative at the origin is given. 

Let r and R < r he the two roots of the polynomial — fxX -|- ft: 

r = fi+ R = gfi — 2/3. 

(See also Q.) From the general theory of differential equations, one has 
Lemma 24. Let v be a solution to 

0 = ^v” -|- /if' -|- /3(f^ — v) (57) 

such that v{x) converges to 1 as x —)• oo. Then, for some non-zero constant A or B, 

• if H one has either 1 — v{x) ~ Ae~^^ or 1 — v(x) ~ Be~^^ as x ^ oo. 

• If gL = -v/2]d = r = R one has either 1 — v{x) ~ Axe~^^ or 1 — v{x) ~ Be~^^ as x ^ oo. 

Furthermore, up to invariance by translation, there are exactly two solutions which converges to 1 
in the fast way (as Be~'^^); one of them approaching 1 by above and the other from below. 
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This lemma simply tells that the solutions to the non-linear equation (57) behave around 
u = 1 as the solutions to the equation linearised around 1. We already used implicitly this result 
in Section [1.2.41 


Proof. This follows from a result of Hartman [12] who shows that if 

X = TX + F{X) 


(58) 


is a non-linear differential system of dimension 2 with T a hyperbolic (eigenvalues have a non-zero 
real part) matrix and F is with F{x) = o(|x|) as x —)• 0 (so 0 is a critical point) , then there 
exists a diffeomorphism (j) with derivative the identity at the origin such that U{t) = (j){X{t)) 
solves 

U = rU. (59) 

Otherwise said the solutions of the linearized system and the solutions of the non-linear system 
are (locally around 0) in one-to-one correspondence through (f). 


We apply this result to the following system where u = 1 — u is a solution of (|57|) 

Xit) = 


x{t) 

y{t) 


u{t) 

u'{t) 


X{t) = 


u'{t) 

u"{t) 


y{t) 

-2ny{t) - 2(3[x{t) - x‘^{t)] 


which has a critical point at {x,y) = (0,0). In this case 


r = 


0 1 
-2/3 -2n 


(60) 


(61) 


with eigenvalues —r = —/i — yf — 2/3 and —R = —y + \/fF — 2/3 (for simplicity we only 


consider the case r ^ R here) and corresponding eigenvectors 
of U = TU are of the form 


1 

—r 


and 


1 

-R 


U{t) = 


Ul{t) 

U2{t) 


= Be 


—rt 


1 

—r 


-F Ae 


—Rt 


1 

-R 


The solutions 


(62) 


Thus the only solutions such that \U{t) \ ~ ce“’'* for some constant c are those such that A = 0. 
If H > 0 then ui approaches by above, if i? < 0 tti approaches by below. Hartman’s Theorem 
tells us that there exists 


^{X)=X + f{X), f{X) = o{\X\)whenx^0 


such that the solutions X{t) of the non-linear system are locally 


xit) = r\uit)). 


(63) 


(64) 


Thus, (after a shift in the argument, replacing x by x -|- In \B\/r) there is exactly one solution X 
to the non linear system such that |X(t)|e’’* has a non degenerate limit and such that x{t), the 
first coordinate of X{t), is eventually positive (resp. eventually negative). □ 


Let s < 1 be such that there is a solution v of (57) with u(0) = s, 1 — v{x) ~ ce“™, and 
v{x) < 1 for all X > 0 (we now know that such an s exists). Then ^^(x) being the maximal 
solution of (Q that starts from s and stays below 1, we must have 

Ws(x) > u(x),Vx > 0. 


Since we also know that the only two possibilities for the asymptotic behavior of uig are that 
either iOs{x)e'^^ —>■ c or ujs{x)e^^ —)> c we conclude that it is the former that holds. The same 
argument apply for ti;s(ic) for any s < sq and in the critical case. This concludes the proof of 
Proposition 1^ for a;(x). 
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2.2.4 Proof of Proposition 


We consider the series <^(2;) = Yln>i defined in (13) with the coefficients an defined in (12). 
The function z 1—is a well dehned object because, by induction on (12) one has easily 
0 < On < 1 and, therefore, TZ > 1. It is then very easy to check by direct substitution that for 
any P G M, the function 


X t-)- v{x) = 1 — ^{Be for x such that |i?|e < TZ, (65) 


is solution to the partial differential equation ]^v" + ^v' + — f) = 0 which appears in (Q 

(when discussing oj) and in ( |26[ ) (when discussing ojg). Recall that r = ji -\- — 2/3 is the 

larger root of + /3. 

As the coefficients an are positive, 4>(2;) is non-negative and increasing for 2; > 0. As ai = 1 
and 02 > 0, it is easy to hnd aO< 2 ;o<l<^ such that <h(zo) > aizo + 02^0 > 1- This implies 
that there must exists a Bo € {0,TZ) (smaller than zq) such that ^{Bq) = 1. With B = Bo, the 
function v{x) in ( |65[ ) is smaller than 1 and converges to 1 for large x as e“™. Using Proposition]^ 
this implies that t)(a:) = ui{x) = 1 — ^{Boe~^^). 

Recall by Theoremj^that Ug for s < 1 is simply equal to uj correctly shifted to have ti;s(0) = s. 
This implies that, for s < 1, iOs{x) = l — ^{Bse~^^) where Bg G (0, Bo] is such that ^{Bg) = 1 — s. 

The case s = 1 is trivial, we now turn to s > 1. As for s = 0, we have the following points: 


for s > 1, Ws is the smallest solution to (26) that remains above 1 (Theorem]^. 


By Lemma 24, there is exactly one solution to (26) which remains above 1 and decays to 
1 as e“™. Because of the previous point, this solution must be uig. 


Now consider <^(2;) for negative arguments. Because <h(0) = 0 and <h'(0) = 1, there must exists 
B G {—TZ, 0) such that <I> is negative on [B, 0). Then, the function x ^ 1 — <h(Re“^*) is solution 
to (26) for s = 1 — ^{B) > 1, remains above 1 for x > 0 and converges to 1 as e“™. Therefore, 
it must be LZg for that particular s. 

But all the functions ujg for 1 < s < so are related through Theorem they are all shifted 
versions of Wsq. Therefore, for any s G (l,so], one has ujg{x) = 1 — ^{Bge~'^^) for a well 
chosen negative Bg (which represents the shift), at least for values of x sufficiently large to have 


\Bg 


< TZ. 


2.3 Proof of Theorem IS 

We assume to be in regime A or B (/x < \/2^) and we want to show how u{t, x) = P^{K{t) = 0) 
converges to a KPP travelling wave. 

The proof is essentially analytic and relies on Bramson’s result [7] and the maximum principle. 
The key step is to compare u{t, x) to a new function : [T, -|-oo) x M 1—)• M (where T > 0 is 
a parameter) where v'^{t,x) is defined as the probability, in the standard branching Brownian 
motion (without absorption nor stopping) with drift /x starting from x, that no particles are 
present in the negative half-line between times t — T and t. In symbols 

v^{t,x) ■.= P^(Vr G [t - r,t],V'U G A/'aii(r) : X„(r) > O), (66) 

where we recall that A/'aii('S) is the population of particles at time s in a branching Brownian 
motion without absorption or stopping. 

The advantage of is that since it is dehned on A^aii it satishes a KPP equation on the 
whole real line: 

dtv'^ = \dxxV^ + + 13{{v^f‘ -v"^), {t,x) G [T,-Foo) x R 

< {T,x) = u{T,x), for X > 0, 

x^(T,x) = 0, for X < 0. 
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Otherwise said the function v^{t,x) = v^{T + t,x) solves the KPP equation on the whole line 
with initial condition i;^(0,x) = u(r, Since for T > 0 hxed, 1 — u{T,x) goes to 0 

as X —)• oo with a super exponential decay (like the tail of a Gaussian), Bramson’s convergence 
Theorem |71 Theorem A], ensures that there exists a constant Ct G K such that we have 


\v^{t, ■ + mt- ^lt + Ct) - 


0 as t —)• oo, 


(67) 


where Bramson’s displacement mt is given in (23). The value Ct depends on T because for 
different T we plug different initial conditions in the KPP equation . 

Since nit — 'nxt-T \/2)dT when t —>■ oo, one obtains taking Ct = Ct — \/2(5T + fiT: 


■ + mt- iJLt + Ct) - 


0 as f —>■ oo. 


Therefore we only need to show that for t large enough, u{t,x) is close to v'^{t,x): 

Lemma 25. 

■) - u{-, •)||oo 0 as T OO. 

In addition, there exists G G M such that Ct —>■ G as T ^ oo. 


( 68 ) 


(69) 


Indeed, assuming that Lemma 25 holds, we can conclude the 


Proof of Theorem\^ Fix e > 0. Using (68) and (69), choose T large enough that ||u^(-, •) — 
u{-, Olloo < e and then choose t large enough so that \\v'^{t, ■ + m-t — fit + Ct) — h*(-)||oo < £■ 
Then, we have that 

||u(t, ■ + m.t- fit) - K{- - G)||oo < \\u{t, ■ + mt- fit) - v^{t, ■ + mt - fit)\\oo 

+ \\v'^{t, ■ + mt- fit) - /i*(- - Gt)||oo 

+ ||/i*(--GT)-h*(--G)||oo, 

< 2e + c|Gr-G|, 

where c = max^-gK h'^{x). As Ct —>• G, for T large enough independently of x this can be made 
smaller than 3e. Thus \\u{t,- + mt — fit) — /i*( - — G)||oo —)• 0 as f —>■ oo, which is the Theorem. □ 


It now remains to prove Lemma We start with 
Lemma 26. For any e > 0 there exists such that for all t > T > one has v'^{t, 0) < e/(1+e). 
(The 1 + e in the denominator makes the following easier.) 


Proof. We use the representation (66). Let t > T; obviously 

v'^{t,Il) < IP’^( min Xu{t) > 0) = h{t,—fit), 

MeA/'all(i) 


where h is the solution of (22). h{t, —fit) is by dehnition the probability that the leftmost particle 
at time t of a driftless branching Brownian motion is to the right of —/it; it is also the probability 
that the leftmost particle at time t of a branching Brownian motion with drift fi is to the right of 
zero. For fi < \/2(d (regimes A and B), this probability is known to tend to zero when t —)• oo. □ 

The next step is the following Lemma: 

Lemma 27. For any e > 0 and any T > one has 

(1 + e)v'^{t, x) — € < u{t, x) < {t, x), {t, x) G [T, oo) x M+. 


(The Tg in Lemma 27 is the same as in Lemma SG) 
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Proof, u < follows immediately from their definitions as probabilities. Let us introduce 

u{t, x) + € 


u{t, x) := 


1 + e 


We have that 


(1 + e)dtu = (1 + e)]^dxxU + (1 + e)ndxU + /3[((1 + e)u - e)^ - ((1 + e)u - e)] 
Performing simple calculations we arrive at 

dtu = ^dxxU + ^idxU + /3{u - l)(u - e + eft) 

> ^dxxU + ^idxU + I3{u - l)ft 

since ft < 1 and e > 0 . 

Now, for any T > T^, we have with Lemma 

0) < ^ ^ = u{t, 0), t>T. 

Moreover one checks directly that 

v^{T, x) = ri(T, x) < u(T, x), x > 0. 

By the parabolic maximum principle (and the unicity of solutions) we get that for any T > 

v^{t, x) < u{t, x), V(t, x) G (T, oo) X M+. 

This proves the first inequality and thus concludes the proof of the lemma. □ 


Lemma 27 implies that \u{t, x) — {t, x)| < e(l — {t, x)) < e for each x G M+ and each t 


and T with t > T > T{€), which is the first assertion of Lemma 25 


The last step is then to prove that Ct has a limit C for large T. 

As u{t, •) is strictly increasing and continuous, u{t, 0) = 0 and lima,_>_|_oo u{t, x) = 1, we may 
define mi : ( 0 , +oo) i—M+ by 

u{t, mi (t)) = 1/2. 

Fix e > 0. We have that 

1 1/2 — h* (mi (t) — mt + fit — Ct) \ < \u{t, mi w) - , mi m 

+ \v'^ (t,mi{t)) — h^{mi{t) — mt + fit — Ct)\, 

< 2 e. 


(70) 


as long as T and t are large enough by ( | 68 [ ) and (69). From this we deduce 
mi{t) -mt + fit-Cr G [/i“^(l /2 - 2 e),/i“^(l /2 + 2 e)]. 

Consequently 

limsup [mi (t) — mt + fit] — liminf [mi (t) — m* + fit] < + 2 e) — hf^{l/2 — 2 e). 

t^+OC 2 t^+OO 2 

Since e can be chosen arbitrarily small we have that limi_^_|_oo [mi (t) — m* + fit] = C, for some 
constant C G M. This and (70) immediately yields that 

lim Ct = C, 

T^ + CXD 


where we used that h^^{l/2) = 0. This concludes the proof of Lemma 25 
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3 Radius of convergence and asymptotic behavior of Sq 


In Section 1.2.3 we related the ujs{x) to a function x i—)• <^( 2 ;) defined as a series of which the 
coefficients an follows the recursive equation (12). We write here the same property in a slightly 
different but equivalent way. Let p € (0,1] be defined by 


p := 


21 

^2 ’ 


and introduce = p^{z/p) and bn^ = Un/p^ These quantities satisfy the relations 

n —1 




1 


n>l 


(n — l)(n — p) 


E'>J "£2- 


(71) 


1 = 1 


Let be the radius of convergence of We know that there exists a i?<jp relating and 
WsQ through 

The following observation will be useful. Since w(,q( 0) = 0 and cn(g(0) < 0, the function 
(defined on a domain containing zero) has a local maximum in zero. This implies that for p > 0 
the function has a local minimum in := pEg^ < 0. In fact is the first local 

minimum (and indeed the first point where the first derivative cancels) one encounters left of 
zero for 

The steps of the proof are the following: 


1. We show that > 4 for small enough p (including p = 0). 

2. We prove that there exists G (—3,0] which is the first minimum one encounters left 
of zero for and that 


['I'^°)]'(x) < 0,x G [a, 771^°^), and ['I'(°)]'(x) > 0, x G (m^°\o], (72) 

for some a G [—3, 

3. We show that converges to uniformly on (—3,0]. This implies that 

lim G (—3,0). (73) 

p\o 


4. Since \pBsq \ —)• < 4 we conclude that that is within the radius of convergence of 

4* for p small enough. The identity 

so = tCso(O) = 1 - ^{Bsq) = 1 - 

shows that 

limpso(p) = 
p\o 

where we made the dependance of sq on p = 2/3/r^ explicit. 

We now prove these points. 

1. The key remark is that if for a real a > 0 and an integer no, one has b^^ < (no — p)a~^ 
for all n G {1,..., no — 1} then, as can be shown by a very simple recursion, the property 
bn^ < (no — p)a~'^ holds for all n > 1. 


25 





Computing the first values of bn \ one checks easily that the maximum of for n G 

14} is around 14.14. For p small enough, by continuity of p ^ bn \ the maximum 
of bn^ for n G {1,..., 14} will be no more than 15 — p and hence one has 

b^n^ < 15 X 4r'^, (for p small enough) (74) 


As a consequence, > 4 for p small enough (including p = 0). 


The bound (74) applies for p = 0. Thus, for any 2 ; G [—3,3] using only the the fifty first 
terms of the expansion leads to an error of at most X]n> 5 i ^ (3/4:)*^ < 3 10“®. In that 
way we computed 3) « —0.8528 and 2.5) ss —0.8575. Therefore 2.5) is 

smaller than both 3) and T(*^)(0) = 0, and the function must have a minimum 

in (—3, 0). In other words we proved G (—3,0). It is easy to check that X]n=i ~ 
l)onX"’“^ > 0.7 for X G [—3,0]. Estimating an error by X]n> 5 i 15n(n — l)(3/4)"' < 0.074 
we conclude that [4'(^Y'(x) > 0 for x G [—3, 0]. In this way we get 


3. By (74) there exist po > 0 and C > 0 such that the functions [T^^^]' are analytic in [—3,0] 


and suppg[o,po], 3 : 6 [- 3 ,o] |'k^^^(x)| < C. By continuity of p 1 —)• bn \ for any x G [—3,0] we have 
[^(p)]'( 

x) —)■ [^(0)]'( 

x). The Vitali-Proter theorem strengthen this to uniform convergence. 
This together with (l7^ implies easily (73). 
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